Anisotropy of magnetothermal conductivity in Sr2Ru04 
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The dependence of in-plane and interplane thermal conductivities of S^RuCU on temperature, 
as well as magnetic field strength and orientation, is reported. We found no notable anisotropy 
in the thermal conductivity for the magnetic field rotation parallel to the conducting plane in the 
whole range of experimental temperatures and fields, except in the vicinity of the upper critical 
field H C 2, where the anisotropy of the H C 2 itself plays a dominant role. This finding imposes strong 
constraints on the possible models of superconductivity in S^RuCU and supports the existence of 
a superconducting gap with a line of nodes running orthogonal to the Fermi surface cylinder. 



PACS numbers: 74.25.Fy, 72.15.Eb, 74.20.-z 

The nature of superconductivity of Sr2Ru04 has at- 
tracted a lot of interest recently |l],^| . An early theoret- 
ical prediction || about the spin-triplet character of the 
superconductivity of the material was supported by sev- 
eral key experimental findings made within the last few 
years. These include (i) the strong dependence of the su- 
perconducting transition temperature T c on the density 
of non-magnetic impurities (ii) the spontaneous ap- 
pearance of magnetic field below T c , as expected for 
a superconductor with broken time reversal symmetry 
(TRS), and (Hi) the lack of any change of the spin sus- 
ceptibility in the magnetic field parallel to the conducting 
plane as expected for a spin-triplet superconduc- 

tor with the Cooper-pair spins parallel to the conducting 
planes. 

To explain experiments (i)-(iii), a p-wave supercon- 
ducting state with broken TRS, described by a d-vector 
d=z(k x + iky), was proposed |J. This state has a full 
gap on the Fermi surface (FS), with the gap mag nitude 
possibly being different on different FS sheets [[l0| . How- 
ever, recent experiments [pd|-p^|, performed on very pure 
samples (characterized by T c close to the estimated value 
for perfect material T c q —1.5 K) revealed a tendency for 
a power-law behavior of quasi-particle excitations down 
to the lowest temperatures (50 mK), contrary to the ex- 
ponential decrease for the fully gapped state. In view of 
these findings it was suggested that d- or /-wave states 
may be better choices for the description of the super- 
conducting state in Sr2Ru04 [jl4|— p~T|| , although the spin- 
singlet d- wave state would be incompatible with the ex- 
perimental result (Hi). 

For a correct description of the superconducting or- 
der parameter in Sr2Ru04, it is important to know the 
location of the line nodes, if there are any, on the FS. 
This information can be obtained only with the use of a 
directional probe. Study of the anisotropy of magnetic 
field dependence of thermal conductivity is one of the 
key experiments for this purpose. Similar experiments 
were performed on the high-T c cuprates, giving a four- 
fold anisotropy consistent with the d-wave order param- 
eter@[l9|. 

In this Letter we report the study of thermal conduc- 



tivity k of Sr2RuC>4 as a function of the orientation of the 
magnetic field within the conducting plane. We studied 
the thermal conductivity with a heat flow direction, Q, 
both along the plane, in [100] direction (ftioo)j and per- 
pendicular to the conducting plane (kooi)- k ioo is elec- 
tronic in origin [^,^l] and thus is directly related to the 
structure of the superconducting gap; however, among 
several factors limiting the accuracy of kioo anisotropy 
measurement, the lack of the symmetry in the magnetic 
field H direction with respect to Q is quite inevitable. 
This makes inter-plane configuration preferable, since 
here H parallel to the planes is always transverse to Q. 
Study of phonon conductivity, which dominates in this 
configuration, offers additionally essential simplification 
of the understanding of quasiparticle density variation. 
In electronic thermal conductivity, generation of quasi- 
particles due to the Volovik effect |2^| should be taken in 
a self-consistent way: on increasing H the increase of the 
density of heat carriers and the decrease of their mean 
free path nearly cancel each other |^||4|. However, the 
density of heat carriers in phonon transport is determined 
only by temperature T and is constant in H, hence all 
the variation with H comes essentially from the increase 
in the density of scatterers, mainly quasiparticles in the 
bulk ||. 

The main result of our study is that we did not find 
any anisotropy beyond the accuracy of our experiment, of 
the order of 0.5%. This is valid in the whole experimental 
range of temperatures T (down to 0.3 K) and magnetic 
fields, except very near to the upper critical field H C 2- A 
two- fold anisotropy up to 2% was found in Kioo, associ- 
ated with the variation of the H direction with respect to 
Q. The upper bounds of the values of the anisotropy are 
essentially lower than theoretically predicted variations 
for the states in which line nodes are running along the 
FS cylinder only, like a d-wave state [ p^|j2rj| . Therefore 
our data are in favor of a state with a line node running 
orthogonal to the FS cylinder, as was proposed recently 

Single crystals of Si'2Ru04 were grown by a Ru self flux 
method in an image furnace p7[ . The samples for our 
study were cut from a surface part of a single crystalline 
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rod, to avoid inclusions of the so-called 3-K phase p8[ . 
The crystals were subjected to annealing to remove resid- 
ual strain caused by cutting. All the samples were charac- 
terized by ac susceptibility and x-ray measurements prior 
to the thermal conductivity study. The data showed the 
high crystalline quality by giving the small width of the 
superconducting transition (10 mK) and the sharp X-ray 
rocking curve peak limited by the instrumental resolu- 
tion. For the measurement of kiooj a sample elongated 
in the [100] direction with the size of 2x0.5x0.05 mm 3 
was prepared. For the Kooi measurement, two samples 
were prepared with the dimensions 0.5x0.5x2 mm 3 . We 
intentionally used samples with the side surfaces running 
in different directions in the conducting plane to check for 
the influence of sample demagnetization factors 29 1. We 
did not find any dependence of the results on the sample 
orientation within the accuracy of our experiment. 

Usual steady state one- heater-two- thermometers 
method was used for n measurements plf| . A miniature 
vacuum cell ]30| ] with the sample was placed in a double 
axis rotator in a 3 He refrigerator insert of a supercon- 
ducting magnet, allowing for a field alignment with an 
accuracy better than 0.1°. The electrical resistivity, p, 
and k were measured with the same electrical/thermal 
contacts. To ensure high resolution of the data it is es- 
sential to keep the thermal gradient within a sample low, 
although this leads to deterioration of the signal to noise 
ratio. To optimize the accuracy under these countervail- 
ing requirements we adopted an intentional overheating 
equal to 0.005 xT, setting the accuracy of our measure- 
ment of k at about 0.5%. We performed measurements in 
both fixed T- H sweep and fixed H-T sweep modes. No 
difference was noticed between these two sets of data. On 
a slow H sweep (typically 1 to 10 mT/min) no hystere- 
sis was observed. High precision H sweep measurements 
take a lot of time and it was difficult to keep the temper- 
ature with the long term stability better than 5 mK in 
varying field and heat load. This limited an accuracy of 
the -ff-sweep data at about 2%. 

In Fig. 1 we show the temperature dependence of 
Kioo/T and kqqi/T. In the normal state, Kioo is pro- 
portional to T, as expected for a normal metal in the T 
range where the scattering of heat carriers is elastic and 
determined by impurities pl[ |. As can be seen from Fig. 
1, the anisotropy of thermal conductivity is much lower 
than the anisotropy of the electrical resistivity. This ob- 
servation indicates the phononic nature of Kooi at 1.5 K. 
In the normal state (in H of 1.5 T parallel to the con- 
ducting plane) the kqqi/T curve is linear down to 0.7 
K, below some curvature appears. The leading k ~ T 2 
term in the temperature dependence of the phonon con- 
ductivity is naturally ascribed to a transport limited by 
scattering on conduction electrons |3l| . 

The electronic contribution to Kooii ^ooij m the normal 
state (in the 1.5 T in-plane field) was estimated in two 
ways. In the first case we compared magnetoresistance in 
pool an d kooi |3§]- Assuming the same relative change of 
both pool and Kqoij this estimation gives Kq Q1 —0. 00947 



Wm 1 K 2 , amounting to 7% of k oi at 1.5 K and 30 
% at 0.3 K. This value almost coincides with nf m de- 
termined from the value of pool by application of the 
Wiedemann- Franz ratio, giving 7% at 1.5 K and 32 % 
at 0.3 K. As we see, in the vicinity of T c , «ooi is mainly 
determined by phonon flow with scattering by conduc- 
tion electrons. In the superconducting state the density 
of quasi-particles decreases and «ooi shows a peak due 
to the phonon mean free path increase. The magnitude 
of the peak, is however, 8 times lower than the expecta- 
tion from the phonon scattering on sample boundaries, 
showing decisive role of scattering by the electronic car- 
riers even in the superconducting state. Simultaneously 
the electronic contribution to k is notably suppressed be- 
low T c pi}] . Therefore, in the superconducting state the 
phonon contribution is dominating, accounting for 88 % 
of kooi at 0.3 K §|. 

Figure 2 shows the field dependence of both kooi and 
Kioo at 0.3 K in a magnetic field aligned along the [100] 
and [110] directions. Despite a clear difference of H C 2 
in the two directions, the slope of Kioo vs. H at H C 2 is 
nearly identical. This observation points to the isotropic 
gap structure near H c i- It is noteworthy that the curves 
for the two heat flow orientations look mirrorlike. This 
fact reflects the opposite effect of quasi-particles on the 
electronic and phononic conductivity, as discussed above. 
Similar curves were obtained with different azimuthal an- 
gles <p of the field within the plane with a 10° step, under 
exact alignment conditions (with the polar angle 0=90° 
adjusted for each measurement). The extracted values 
of the kooi at several representative fields, as well as the 
anisotropy of the upper critical field, are shown in the 
polar plot of Fig. 3. As can be seen, H C 2 shows a no- 
table four-fold anisotropy, with the maximum and the 
minimum in the [110] and [100] directions, respectively. 
The anisotropy of H C 2 amounts to 4 %, in a reasonable 
agreement with the AC susceptibility (3^] and the in- 
plane thermal conductivity pl[ ] data. The scatter of «ooi 
values in this measurement is about 2%. 

The difference of the magnetic field effect should be the 
largest for the two principal directions, [100] and [110], 
within the conducting plane. In view of the tetragonal 
symmetry of the crystal the superconducting gap can 
have nodes only in one of these directions. As can be 
seen from the comparison of the curves in Fig. 2, the dif- 
ference, if any, is very small. To eliminate the problem of 
temperature drift during H sweep (which is the essential 
noise factor) we concentrated on measurements in the T 
sweep mode. We performed the measurements with fixed 
H from 0.1 to 1.2 T, with a 0.1 T step, for the [100] and 
[110] directions. It is essential to cover such a broad field 
range, since different effects can affect the anisotropy at 
low and high fields. In Fig. 4 we show the results of these 
measurements at some representative fields, presented as 
deviations from the average k value. We find that the 
deviations are random and less than 0.5 %. 

In order to understand this result, first we recall the 
anisotropy in the quasiparticle density generated with 
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field, predicted for a d-wave state by Vekhter et al. [ p4| , 
with the difference for [100] and [110] directions amount- 
ing up to 30 %. To our knowledge there is no experi- 
mental observation of this effect. For the cuprates the 
limitation in experimentally accessible values of the re- 
duced magnetic field, H/H c2 <S 1 is one of the possible 
reasons. For Sr2Ru04, in which the whole field range 
is accessible, the anisotropy of the H C 2 itself is making 
a dominant contribution in the very vicinity of H c i- At 
1.2 T, near (but not very close) to H c2 , the anisotropy 
is already very small. Since the interplane heat flow is 
influenced only by the quasiparticle density, the latter 
remains constant within the accuracy of our experiment. 
Thus such anisotropy effect is not present here. Another 
contribution is the anisotropy in the scattering of quasi- 
particles due to the Andreev reflection on the superfluid 
current around the vortex core. This effect would be 
important for electronic contribution at low fields and 
is responsible for the weak-field anisotropy observed in 
the cuprates p8| , ^9| . Again, there is no evidence for this 
effect in Sr 2 Ru04. 

In this Letter, we evaluated the anisotropy by covering 
the whole range of the field. The lack of the anisotropy 
both at low and high fields makes it reasonable to con- 
clude that both of the abovementioned mechanisms are 
not operating here. The isotropic quasi-particle contri- 
bution can come either from a line node in the super- 
conducting gap running orthogonal to the FS cylinder, 
or from a sheet of the FS whose superconducting gap is 
easily and isotropically destroyed by the in-plane field. 
The former situation is quite analogous to the case of 
UPt3, in which the presence of the horizontal node is a 
likely reason for the lack of anisotropy |3^] . The latter is 
a natural extension of the orbital - dependent supercon- 
ductivity (ODS) model |l(]] . From the anisotropy exper- 
iment alone it is not possible to make the preference for 
either the ODS or the "horizontal" node model. However, 
for a full gap state it is difficult to explain the power-law 
dependence of the quasiparticle density at low tempera- 
tures |lj,[ll|. Therefore, an ODS model with full gaps on 
all three FS sheets does not give an appropriate descrip- 
tion of the complete set of experiments. Measurements 
at lower temperatures may clarify this issue. 

In conclusion, we have found no notable anisotropy of 
the thermal conductivity in Sr2Ru04 with high experi- 
mental resolution. In conjunction with the existence of 
low temperature quasi-particle excitations, indicative of 
nodes in the gap, our results strongly suggest the pres- 
ence of the line node running orthogonal to the Fermi 
surface cylinder. 

The authors acknowledge A. E. Kovalev and V. A. 
Bondarenko for the help in experiment design, and D. 
Agterberg, A. G. Lebed, K. Maki and M.Sigrist for valu- 
able discussions. 

Note added. In the following Letter, K.Izawa et al. 
p6| make an important extension of our previous study 
of anisotropy in kioo plfl • 
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FIG. 1. Temperature dependence of the in-plan 
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FIG. 2. Field dependence of thermal conductivities Kooi 
and ftioo at 0.3 K for the magnetic field along two principal 
directions, [100] and [110], in the conducting plane. Insets 
show an expanded view of a region near H rmc 2 ■ 



4 



